Cosmic Bose dark matter 
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Cold and hot dark matter (CDM, HDM) imprint distinctive effects on the cosmological observ- 
ables, naturally, they are often thought to be made of different kinds of particles. However, we 
point out that CDM and HDM could share a common origin, a Bose-Einstein condensate emerging 
with two components. In this framework, the mass and temperature constraints on HDM contain 
fundamental information of CDM, but also, the critical temperature of condensation shall be larger 
than the HDM temperature. We discuss two scenarios: a gas made of bosons and a gas made of 
boson-antiboson pairs. We use some cosmological data surveys to test the idea and constrain the 
bosonic DM parameters, including a forecast for the Planck mission. We find that the bosonic DM 
picture is consistent with data in the first scenario, although, the second one might be increasingly 
interesting for future data surveys. 
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One of the most intriguing questions in modern cos- 
mology comprises the nature of the so far unidentified 
one quarter part of the Universe content, known as dark 
matter (DM) . Cosmological observations seem to support 
the idea that DM is made of some kind of non baryonic, 
non relativistic, and weakly interacting massive particles. 
Many efforts to elucidate this issue have been taken over 
the last decade, mainly motivated by the idea that the 
answer to the riddle will very likely change our current 
understanding of matter and its interactions. One fa- 
vorite alternative is that DM could be sourced by a Bose- 
Einstein (BE) condensate made of light particles. Several 
noteworthy candidates have been addressed in this cate- 
gory e.g., axions and axion-like particles scalar fields 
[2] and perhaps more exotic species such as bosonic neu- 
trinos and massive photons [3| . The landscape of bosonic 
DM appears attractive enough for a deliberate study and 
test of its extents of validity. Our pursuit is not to an- 
alyze one specific candidate, but instead, outline generic 
properties of bosonic DM and explore their constraints 
with analysis of cosmological data. 

At very early stages of the Universe, DM is thought to 
be in thermal contact, either with itself or other species, 
composing a primeval fireball. DM candidates can be 
classified according to their particle velocities dispersion, 
for our purposes: (1) Particles with vanishing velocity 
dispersion are termed cold dark matter (CDM). (2) Hot 
dark matter (HDM) are particles that, being ultrarela- 
tivistic at early stages, become nonrelativistic only at re- 
cent epochs. The standard candidate for HDM is the neu- 
trino, whose decoupling thermal energy is at the MeV's 



scale and the temperature today is around 1.95 K. Dur- 
ing the expansion of the Universe, while HDM is free- 
streaming, CDM is gravitationally clustering. A Uni- 
verse dominated by HDM would predict a suppression 
of small matter structures, in contradiction with mea- 
surements. Hence, CDM is favored to form the overall 
structures and compose most of the DM in the Universe. 
These two families are usually assumed to be unrelated 
to each other; nevertheless, the DM paradigm docs not 
preclude scenarios that involve jointly CDM and HDM. 
Let us now introduce our assumptions framework, and 
call it the bosonic DM picture: DM is dominated by par- 
ticles in a BE condensate state. Being relativistic at de- 
coupling from the primeval fireball, the ground state is 
composed of condensed bosons which act as CDM, but a 
complementary fraction of thermal bosons resides in ex- 
cited states. As the Universe expands the temperature 
decreases, then, thermal bosons might pass through a 
nonrelativistic transition, being present today as bosonic 
HDM. Following this idea, we describe the full bosonic 
system in equilibrium as the sum of two momenta dis- 
tributions: thermal bosons obey the standard BE mo- 
menta distribution, while that of condensed bosons ex- 
hibits a delta-like shape. This allows us to evolve ther- 
mal bosons in phase-space and condensed bosons as a 
fluid. Bosonic DM evolves along with baryons, photons, 
and other species, coupled only gravitationally. These 
dynamics determine the observed matter structures and 
the cosmic microwave background (CMB) pattern, which 
can be tested with cosmological data today. We discuss 
two scenarios likely to yield relativistic BE condensation: 
a gas formed purely of bosons Q and a gas of boson- 
antiboson pairs We analyze both through their ex- 
pressions for the critical temperature of condensation. 
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thereby, testing their extents of consistency as DM can- 
didates. Given the distinctive effect of CDM and HDM, 
the bosonic DM picture might be a way to infer infor- 
mation of CDM, hke the mass, extracted from thermal 
bosons, but also, their temperature should be related to 
interaction couplings of CDM particles. 

Let us begin to outline the BE condensation inside 
the primeval fireball, at 2:ero-order in perturbation the- 
ory. The process of condensation should be driven by 
self interactions of high-energy bosonic particles Q . Ac- 
cordingly, bosons obey the BE momenta distribution 
/^°^(p) g^/ (e(^"^)/^* - l) in the relativistic limit 
TO << T^, where to is the mass of bosons, their tem- 
perature, fj, the chemical potential, and is the num- 
ber of internal degrees of freedom, /be defines the crit- 
ical temperature of condensation Tc- In the relativistic 
regime, if the bosons, 4>, are their own antiparticle, then 

T^^ = (7r2n(°)/g0C(3))^^^ 0, where C(3) w_1.202; but, 
in the case of abundant boson-atiboson (fx/) pairs pro- 
duction T;!'^ = (3An(o)/50TO)^^^[5l. Here, n^") is the 
total number density, and An^^^ the excess of bosons 
over antibosons number density. In either case, at tem- 
peratures below the critical temperature and the limit 
/i — >■ TO, the theory predicts a phase transition in which 
a large fraction becomes allocated in the lowest energy 
occupation numbers, forming a coherent state. Hence- 
forth, the momenta distribution hlows-up, taking a singu- 
lar shape; when thermal equilibrium is reached, it reads 
/(o)(p) = ni°^(53(p) + /^°j^(p) 0, where is the number 
density of condensed bosons. 

In the instantaneous decoupling approximation, the 
subsequent evolution of the full bosonic system is af- 
fected only by the expansion rate of the Universe and 
small gravitational instabilities. These instabilities are 
expressed as linear perturbations to the space-time ge- 
ometry. In conformal newtonian gauge, the metric reads 
ds2 = a(r)2 [-(1 + 2iP)dT'^ + (1 - 2ip)dxUxi] . Here a 
is the scale factor of the Universe's expansion, r the 
proper time, and ip along with ip are scalar modes. We 
describe the effects on the dynamics of thermal bosons 
by the usual approach to HDM 0]. The momenta dis- 
tribution is perturbed to first order as /t(x, p,r) = 
/(0)(p)[1 + v[/(x,p,t)], with - 5./(e^/^* - !)■ 

Then, is expanded in a Legendre series to obtain a 
hierarchical system of coupled Vlasov equations. On the 
other hand, the phase spaece distribution of condensed 
bosons may be written as a narrow Maxwell-like one Q, 



/,(p)== (2^)3/2^ exp 



(p + mv) 



(1) 



where v = |v| is the linear perturbation bulk flow speed. 
The width a of the distribution is certainly small, mani- 
festing the distinctive blow-up of BE condensation. The 
spatial overdensities of condensed bosons are dcflned in 



for /c, two fluid equations can be obtained for S and 
V. In Fourier space they read 6k + ikvk ~ 3ipk, and 
Vk + = ikipk [lo| . where k = |k| is the magnitude 
of the Fourier wave-number vector. These are formally 
the dynamic equations of CDM, often obtained assum- 
ing very massive particles. In contrast, condensed bosons 
are expected to be very light; but also, the narrow shape 
of fc makes negligible all quadratic terms of p/E ~ p/m 
for the condensed fraction. Hence, under this framework, 
even light particles could behave as CDM, provided they 
posses small velocities dispersion, which is attainable by 
the coherent state of condensed bosons. 

The evolution of bosonic DM after decoupling is as fol- 
lows. Their velocity and temperature can be affected only 
by the expansion rate of the Universe. Also, in a refer- 
ence frame, comoving with the expansion of the Universe, 
the total number density of bosons is conserved. Conse- 
quently, the number density of particles evolves like 
and the temperature, as well as the critical temperature, 
like . If thermal bosons become nonrelativistic at 
recent epochs (to > T^), the zero-order energy density 



reads pnr 



TO ), where nx is the number den- 



sity of thermal bosons and depends on their temperature. 
Thus, the content of bosons has two nonrelativistic com- 
ponents today: the content of bosonic CDM (B-CDM) 
ilc and bosonic HDM (B-HDM) fJ//, which can be si- 
multaneously estimated using cosmological data. This 
is an interesting scenario because cosmological data can 
provide constraints on the mass of bosonic DM by means 
of 
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terms of Ur 
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(1 -I- 6). From the Vlasov equation 



where pcr is the critical density of the Universe as mea- 
sured today in terms of the Hubble constant Hq [lH . 

In Eq. ([2]), we need the temperature in order to know 
the mass of bosonic DM. We can obtain it by recalling 
that the Universe expands adiabatically, i.e., the entropy 
remains constant in a comoving reference frame. At very 
early epochs, when the Universe is dominated by radia- 
tion, the entropy is s = (27r2/45) g^{T)T^ [H], where T 
is the temperature of the primeval fireball, and (7* is an 
effective number of internal degrees of freedom of all the 
relativistic species in thermal contact. (CMB photons 
are the last species to decouple, thus, today g^xi = 2, 
corresponding to two photon polarization states.) The 
entropy leading order term of the bosonic gas is ^ T^, 
as a consequence, the temperature of thermal bosons is 
proportional to that of CMB photons. This proportion 
depends on the number of relativistic degrees of freedom 
of the thermal bath from which they decouple, which 
could include standard and even some exotic species. To 
ilustrate, after neutrino decoupling (T < 1 MeV), elec- 
trons and positrons start to annihilate into photons, 
leading to an increase in the temperature of photons with 
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respect to neutrinos by a factor proportional to the inter- 
nal degrees of freedom of electrons and positrons. Conse- 
quently, the temperature of neutrinos today (and thereby, 
the usual temperature of HDM) is T^^ = (4/11)^^^ Tcmb, 
where Tcmb is the measured temperature of the CMB. 
For bosons, we model any deviation from the usual HDM 
temperature by means of 



T, 



4>,o 



11 

T 



.9x 



-1/3 



^cinb ; 



(3) 



where gx is negative (positive) if bosons decouple during 
(before) the process of annihilations. A particular 
model of bosonic DM in this scheme would contain a 
relation of with their interaction couplings, here, we 
regard g^ just as a free parameter. An increase of 5x 
reduces the free-streaming of bosonic HDM, yielding less 
suppression of small scale matter structures, resembling 
the effect of CDM. 

As an additional remark, we address a bound on the 
strength of bosons self interactions. It comes from the 
fact that our study relies on the range of validity of 
the dilutcncss approximation Q, which can be written 
for condensed bosons in terms of the scattering cross 
section, ascatt ^ (ni"'')^^/^. Besides, as a minimal 
condition, the content of B-CDM should be at least 
larger than baryonic; we find ascatt ^ 87r(m/pcrf^&)^^'^, 
which is more restrictive for ultralight models. Con- 
versely, for B-HDM it may be written as Ugcatt ^ 
3/4 + gx)^H I {g^^b)f^ cm\ 

In the generic study we are presenting here, the bosonic 
DM parameters are addressed as the mass, m, and the 
factor bounds on their values can be obtained from 
a statistical analysis on cosmological data. To this end, 
we perfom a simultaneous analysis [l^l to data pro- 
vided by the 7-year WMAP, ACBAR and QUaD surveys 
for the CMB temperature anisotropics 11411. together with 
the SDSS large redshift galaxy survey [15[. We will call 
this set W-I-AQS. The input model fit to data is obtained 
from an adapted version of the public code CAMB (lij . 
which includes our bosonic DM component. We study 
two models of scalars {g^ = 1) given by the case of a gas 
of bosons and boson-antiboson (jxj) pairs; in the rela- 
tivistic limit, they differ only by a factor of two in the 
momenta distribution of thermal bosons. For simplicity, 
we consider a flat geometry with cosmological constant 
and choose to vary the following set of cosmological pa- 
rameters: the content of baryons Jib, the content of CDM 
He, the scalar inital amplitude A^, the spectral index n^, 
the optical depth of reionization , and one free bias pa- 
rameter for the SDSS data. We make use of flat priors 
with starting points at g^ = 0, m = 0, and the best- 



fit values found by the WMAP team [17[ for the rest of 
parameters. 

The posterior probabilities of bosonic DM parameters 
are ploted in Fig. [T] and listed in table HI from which 
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FIG. 1: 1,2, 3-a probability contours from WMAP-I-AQS for 
the bosonic HDM parameters in the case of bosons (up) and 
boson-antiboson pairs (down) models. Blue contours refer to 
Planck-hAQS. 



follows constraints in the temperature of bosons Tg = 
2.14±0.02 K, and boson-antiboson gas T^^ = 1.91±0.05 
K (95%). The bound in the cross section of self interac- 
tions becomes ascatt ^ 0.18 cm^. From Eq. ([2]), we may 
set lower limits on the critical temperature of condensa- 
tion (evaluated today): T^^q > 9.53 K and T^^^ > 1.91 
K, for the two cases. It is interesting to explore the con- 
straints that the Planck mission will achieve on our model 
parameters. We take a fiducial power amplitude given by 
the best fit to the WMAP data, we neglect systematic 
effects and assume a sky coverage of 65%. The data fore- 
cast is based on the contributions from cosmic variance 
and noise specifications of the first three channels of the 



(95 %) 


3x 


m (eV) flc X 10"' 


(l> W+AQS 


-1.42±0.11 


<0.14 2.1±0.025 


P-I-AQS 


-1.42±0.05 


<0.11 2.1±0.010 


<j)<j) W-I-AQS 


0.67±0.26 


<0.26 2.5±0.030 


P-f AQS 


0.67±0.16 


<0.16 2.5±0.012 



TABLE I: Posterior probabihties of bosonic DM parameters 
from the two data sets considered, for the case of bosons (j> 
and boson-antiboson dxh. 
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FIG. 2: Critical temperatures of condensation (evaluated to- 
day) as a function of the mass for each model. For bosons, 
the lower limit T^q > 10.33 K is higher than our best fit of 
Tq = 2.14 K. For boson-antiboson, the limiting-critical tem- 
perature is comparable with the bosonic temperature. 



bounds on the critical temperature of a boson-antiboson 
gas. Our forecast for Planck data predicts important im- 
provements in the parameters estimation. Additionally, 
we observe that thermal bosons and standard neutrinos 
are both compatible with the presented data. This sys- 
tematic along with the properties of bosonic DM could 
be further studied with future and larger cosmological 
surveys. 
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dence. IRM thanks the kind hospitality of the APC, 
Paris Diderot University, where a stage of this work was 
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high frequency instrument (see |18|); we call this set of 
data as P+AQS. As appreciated in Fig. [1] and table U 
our Planck forecast reduces considerably the probability 
parameter space. Subsequently, we find Tq = 2.14±0.009 
K, >10.33 K for bosons and T;}^ ^ 1.91±0.009 K, 
T^^Q > 1.92 K for boson-antibon, the bound on the self in- 
teraction cross section becomes ascatt ^ 0.15 cm^ (95%). 
Our result Tq < 0.22 T^q indicates consistency with the 
assumptions of the bosonic DM picture. Even if the case 
of boson-antiboson is not conclusive in this respect, note 
however, that the critical temperature depends mainly 
on the mass constraints. We show this dependence in 
Fig. [5] with the temperature fixed to our best fit. A 
boson-antiboson gas exhibits a huge critical temperature 
if their particles are ultralight. Therefore, better con- 
straints on the mass of bosons by larger surveys could 
yield upper bounds on the critical temperature, perhaps 
in consistency with the hypothesis of BE condensation. 

Finally, from a similar analysis for fermionic HDM (F- 
HDM) we find bounds in the sum of neutrino masses and 
the number of extra relativistic species, ^ nii, <0.45 eV, 
iV=1.10±0.18 (95%), in concordance with some previous 
reports ^ . When we compare the F-HDM and B-HDM 
spectra, differences are roughly apparent at the ten per- 
cent level; however, our x^-study shows no signficant dif- 
ference between B-HDM and F-HDM in their ability to 
fit the presented data. This should be considered as a 
source of systematic uncertainty in cosmic neutrino pa- 
rameter estimations [l9| . 

In summary, we have presented a generic study of 
DM based on BE condensation, from which, CDM and 
HDM, are intrinsically related. From our statistical anal- 
ysis, we constrain the mass of bosonic DM (which in our 
framework is the same for both CDM and HDM), the 
temperature, and the self interaction cross-section. We 
find that the bosonic DM picture is compatible with the 
data used here for a gas of bosons, and discuss how fu- 
ture constraints on the mass would yield much higher 
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